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ABSTRACT 


The  Free  Electron  Laser  has  the  potential  to  become  a  revolutionary  weapon 
system.  Deep  magazines,  low  cost-per-shot,  pinpoint  accuracy,  and  speed  of  light 
delivery  give  this  developing  weapon  system  significant  advantages  over  conventional 
systems.  One  limiting  factor  in  high  energy  laser  implementation  is  thermal  blooming,  a 
lensing  effect  which  is  caused  by  the  quick  heating  of  the  atmosphere,  so  that  the  laser 
beam  does  not  focus  on  the  desired  spot,  thereby  degrading  the  effectiveness  of  the  laser 
on  target.  The  use  of  multiple  beam  directors  focusing  on  a  target  from  a  single  platform 
may  mitigate  thermal  blooming  by  allowing  half  of  the  laser’s  energy  to  travel  through  a 
given  volume  of  air,  so  that  they  only  overlap  very  near  the  target.  Less  energy  traveling 
through  a  given  volume  of  space  means  less  heating,  and  therefore  lessens  the  effects  of 
thermal  blooming.  Also,  simulations  of  LEL’s  were  conducted  modifying  parameters 
such  as  the  number  of  undulator  periods,  electron  beam  focus,  the  normalized  Rayleigh 
length,  and  mirror  output  coupling,  in  order  to  determine  optimum  design  parameters. 
New  parameters  for  the  next  proposed  EEL  were  simulated  to  examine  the  effect  of 
mirror  tilt  on  laser  power  and  extraction  as  well. 
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I.  INTRODUCTION 


A,  JUSTIFICATION  FOR  THE  FREE  ELECTRON  LASER  (FEL) 

Over  the  history  of  naval  warfare,  there  have  been  eertain  teehnologieal 
developments  that  have  been  touted  as  revolutions  in  military  affairs.  Among  these  ean 
be  ineluded  the  use  of  the  sail,  eanons,  steam  power,  metal  hulls,  and  the  evolution  of 
missile  and  radar  teehnology.  The  next  step  in  the  evolution  of  naval  warfare  may  be  the 
use  of  eleetrie  weapons,  sueh  as  lasers  and  rail  guns.  A  strong  eontender  for  naval 
applieation  of  direeted  energy  is  the  free  eleetron  laser.  It  has  eountless  applieations  and 
advantages  over  other  laser  teehnologies,  whieh  will  be  discussed  in  this  section. 

1,  Laser  Applications  in  the  Navy 

Lasers  have  many  applications  in  the  maritime  arena.  The  first  and  most 
prominent  is  their  use  in  anti-ship  cruise  missile  defense.  Self-defense  weapons,  such  as 
the  CIWS  (Close-In  Weapon  System)  and  RAM  (Rolling  Airframe  Missile),  currently  in 
the  U.S.  arsenal  are  insufficient.  Both  have  very  limited  range,  difficulties  intercepting 
maneuvering  missiles,  and  limited  magazines.  Lasers  possess  the  potential  to  have  much 
longer  ranges,  virtually  bottomless  magazines,  and  fewer  complications  when  attacking  a 
maneuvering  target. 

The  laser’s  range  is  limited  by  the  line  of  sight,  making  ranges  of  8  to  10  nm 
possible  in  good  weather.  The  CIWS  has  a  range  of  less  than  2  nm  and  the  RAM  6  nm. 
The  laser  adds  to  this  range  advantage  the  capability  of  zero  flight  time.  Because  the 
laser  propagates  at  the  speed  of  light,  the  effect  begins  instantaneously,  whereas  for  the 
other  systems,  there  is  a  significant  delay  between  launch  and  intercept.  This  delay 
between  engagement  and  kill  can  be  deadly,  as  it  leaves  little  reaction  time  to  fire 
additional  salvos  if  the  first  is  ineffective. 

The  term  “bottomless  magazine”  is  applied  to  electric  lasers  because  they  do  not 
use  conventional  ammunition.  Electric  laser  technologies  are  limited  only  by  the 
electricity  sent  to  it,  which  depends  on  the  ship’s  electrical  power  available.  As  long  as 
the  ship  has  the  fuel  to  supply  electric  power  to  the  laser,  the  laser  operates. 
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Lasers  are  expeeted  to  be  mueh  more  effective  against  high-speed  targets  capable 
of  high-G  turns.  This  is  because  there  is  no  longer  the  need  to  lead  the  target.  With 
current  conventional  weapon  systems,  there  is  the  aforementioned  time  lag  between  the 
time  a  round  is  fired  and  when  it  makes  contact  with  the  target.  As  such,  the  weapon 
system  must  lead  the  target,  shooting  where  it  predicts  the  target  to  be  once  the  rounds 
meet  the  missile.  If,  however,  the  missile  makes  too  many  maneuvers  during  this 
intervening  time,  the  rounds  can  miss.  Similarly,  when  using  missiles,  the  missile  must 
attempt  to  hit  a  target  missile  moving  at  an  incredibly  high  relative  speed.  If  the  target 
missile  maneuvers,  it  is  difficult  for  the  missile  tracking  it  to  maneuver  quickly  enough  to 
destroy  it.  The  advantage  of  speed-of-light-delivery  is  that  all  of  these  considerations  are 
now  insignificant.  The  laser  director  needs  only  to  adjust  the  tilt  of  a  mirror  to  maintain 
track,  and  does  not  need  to  attempt  to  lead  the  target.  No  matter  what  the  target  missile 
does,  it  cannot  out-maneuver  the  laser. 

There  are  numerous  other  applications  for  lasers  in  the  military.  One  of  these  is  a 
capability  which  the  current  arsenal  is  lacking,  the  “soft  kill”.  A  soft  kill  is  when  a 
system  is  rendered  ineffective  due  to  means  other  than  the  actual  structural  damage 
usually  associated  with  weapon  systems.  The  laser  can  accomplish  many  variations  of 
the  soft  kill.  One  is  sensor  blinding,  in  which  the  target’s  means  of  detection  and  tracking 
are  disabled.  The  laser  can  accomplish  this  by  overloading  optical  circuits  or  destroying 
key  electronic  components.  Another  is  glazing,  in  which  the  laser  is  aimed  at  a  glass 
surface  and  the  interaction  causes  the  glass  to  cloud  as  in  Figure  1  below.  This  can  be 
extremely  effective  against  aircraft  cockpits.  A  third  method  of  soft  kill  is  the  ability  to 
disable  vehicle  engines.  Lasers  can  pin-point  specific  areas  of  vehicles  and  burn  through 
the  body  in  order  to  destroy  key  engine  components,  especially  electrical  components 
such  as  spark  plug  wires. 
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Figure  1.  Canopy  Glazing  by  Laser 

This  capability  leads  to  another  aspect  of  the  laser  that  far  surpasses  any  weapon 
eurrently  in  the  inventory,  minimum  collateral  damage.  In  today’s  eonflicts,  there  is  an 
emphasis  on  aceuracy.  This  has  led  to  the  implementation  of  precision-guided  bombs  as 
the  munitions  of  choice.  As  aeeurate  as  is  eurrent  laser  and  GPS  guiding  teehnology, 
there  is  still  a  delicate  balance  that  must  be  maintained.  An  explosive  should  be 
minimized  and  the  blast  aimed  in  order  to  minimize  eollateral  damage,  but  still  large 
enough  to  produce  the  desired  lethality.  There  is  a  limit  on  how  precise  an  explosive  can 
be  due  to  the  blast  effects.  With  lasers,  the  energy  is  localized  on  a  small  area  of  the 
target,  so  that  there  is  little  effect  outside  of  the  aim  point.  This  allows  for  unprecedented 
accuracy  and  minimal  damage  outside  of  the  desired  area. 

Finally,  there  are  applications  for  the  laser  in  the  field  of  ship’s  self  defense.  This 
applies  to  such  scenarios  as  swarm  raids,  in  which  an  adversary  uses  a  large  number  of 
small  boats  such  as  speed  boats  or  jet  skis  armed  with  hand  held  rockets  and  small  arms 
in  order  to  overwhelm  the  ship’s  defenses.  Currently,  ships  are  equipped  with  a  few  .50 
caliber  machine  guns  and  perhaps  a  20mm  chain  gun.  Cruisers  and  destroyers  also  have 
their  5”  gun  mounts.  Fire  hoses  are  used  to  repel  boarders.  All  of  these  systems  are 
inaeeurate  and  diffieult  to  use  against  such  high-speed  threats,  and  there  is  slow  transition 
from  one  target  to  the  next.  A  laser  would  not  be  affeeted  by  the  high-speed  maneuvers 
of  which  these  platforms  are  capable  for  the  same  reasons  that  it  is  not  affected  by 
missiles’  maneuverability.  The  laser’s  accuracy,  speed  of  light  delivery,  and  ability  to 
quickly  shift  from  target  to  target  are  immeasurably  better  than  those  currently  in  use. 


3 


This  aspect  of  the  laser’s  potential  use,  as  well  as  the  minimum  collateral  damage 
and  soft  kill  capability,  can  also  be  applied  to  unknown  threats  approaching  a  ship. 
Currently,  the  ship’s  force  has  few  options  other  than  lethal  force  to  defend  itself  In  the 
current  world  climate,  threats  are  not  usually  well  defined.  As  evidenced  by  the  USS 
Cole  incident,  there  is  little  a  ship  can  do  but  wait  to  find  out  the  intentions  of  an 
approaching  unidentified  vessel.  If  lethal  force  were  used  to  deter  an  approaching  vessel 
that  had  not  identified  itself,  the  ramifications  would  be  severe  if  it  was  a  civilian.  Lasers 
would  allow  ship’s  force  to  disable  an  approaching  vessel  whose  threat  status  could  not 
be  determined. 

2,  Advantages  of  the  FEL  Over  Other  Laser  Technologies 

The  FEL  has  many  advantages  over  the  other  laser  technologies,  solid-state  and 
chemical,  available  today.  These  include  scalability,  the  property  of  being  tunable,  and 
the  lack  of  hazardous  chemicals  and  materials  as  output  products. 

Free  Electron  Easers  are  currently  the  most  scalable  lasers  available.  The 
potential  for  mega-watt  operation  is  much  more  feasible  for  FEE’s  than  any  other 
medium.  This  is  primarily  because  the  lasing  medium  is  a  vacuum.  In  all  other  lasers, 
there  is  a  medium  in  which  lasing  occurs  that  must  be  cooled,  and  this  puts  a  significant 
constraint  on  the  amount  of  power  that  can  be  achieved.  In  solid  state  lasers,  there  is  a 
solid  crystalline  medium  through  which  the  laser  beam  must  propagate  in  order  to 
increase  its  power.  Absorption  can  significantly  heat  the  material,  so  that  it  must  be 
cooled.  Because  of  the  solid  nature  of  this  material,  the  heat  must  be  conducted  out, 
which  is  limited  by  the  surface  area  that  can  be  reached  by  the  cooling  medium. 

With  a  chemical  laser,  this  lasing  medium  is  exhausted,  and  the  heat  carried  out 
with  this  exhaust.  While  this  seems  to  mitigate  the  cooling  problem,  this  exhaust  gas  is 
highly  toxic  and  corrosive.  In  a  ship-board  environment,  this  exhaust  is  very  difficult  to 
direct  such  that  it  is  not  detrimental  to  personnel  and  equipment.  Also,  the  storage  of 
these  lasing  gases  is  a  serious  issue,  since  the  chemicals  are  toxic  and  corrosive  before 
reacting  as  well. 

A  tunable  FEE  means  that  it  has  the  capability  of  changing  the  wavelength  of  the 
light  emitted,  or  tuning,  as  required,  without  significantly  impacting  the  design  of  the 
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laser.  Solid  state  and  chemical  lasers  depend  upon  the  materials  used  to  determine  their 
operational  wavelength.  In  order  for  solid  state  or  chemical  lasers  to  laze  at  different 
wavelengths,  new  mediums  must  be  found  with  the  appropriate  properties,  and  the  whole 
laser  effectively  redesigned.  Tunability  is  useful  because  the  atmosphere  absorbs  some 
wavelengths  of  electromagnetic  radiation,  while  allowing  others  to  pass  relatively 
unaffected.  Due  to  its  ability  to  adjust  the  wavelength  with  little  impact  on  design,  PEL 
designers  can  look  for  what  wavelengths  best  propagate  through  the  atmosphere,  and 
build  the  appropriate  design  to  best  take  advantage  of  these  windows. 

B,  THESIS  CONTENTS 

This  thesis  will  begin  by  discussing  the  theory  behind  a  free  electron  laser.  It  will 
then  discuss  a  possible  solution  to  thermal  blooming,  which  is  a  prominent  limiting  factor 
in  high  energy  laser  implementation.  Thermal  blooming  is  a  lensing  effect  which  is 
caused  by  the  quick  heating  of  the  atmosphere,  so  that  the  laser  beam  does  not  focus  on 
the  desired  spot,  thereby  degrading  the  effectiveness  of  the  laser  on  target.  This  potential 
solution  is  the  use  of  multiple  beam  directors  from  a  single  platform  to  focus  on  a  target. 
This  may  mitigate  thermal  blooming  by  allowing  half  of  the  laser’s  energy  to  travel 
through  a  given  volume  of  air,  so  that  they  only  overlap  very  near  the  target.  Less  energy 
traveling  through  a  given  volume  of  space  means  less  heating,  and  therefore  lessens  the 
effects  of  thermal  blooming. 

Also,  simulations  of  FEL’s  were  conducted  modifying  parameters  such  as  the 
number  of  undulator  periods,  electron  beam  focus,  the  normalized  Rayleigh  length,  and 
mirror  output  coupling,  in  order  to  determine  optimum  design  parameters.  Simulations 
were  also  conducted  in  order  to  examine  the  effect  of  mirror  tilt  on  laser  power  and 
extraction. 
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II.  THEORY 


A.  BASIC  FEL  DESCRIPTION 

The  Free  Electron  Laser  may  be  described  in  completely  classical,  as  opposed  to 
quantum-mechanical,  terms.  It  is  easy  to  explain  qualitatively,  but  the  quantitative 
physical  description  is  challenging  enough  to  be  interesting.  The  FEL  is  based  on  the  fact 
that  accelerating  electrons  radiate  energy.  In  order  to  achieve  this,  electrons  are  emitted 
from  a  cathode  in  discrete  bunches.  These  are  accelerated  in  an  RE  accelerator  to 
relativistic  velocities.  These  electrons  are  then  sent  into  what  it  called  an  undulator, 
which  contains  a  stationary  magnetic  field  in  a  periodic  transverse,  spatial  configuration. 
There  are  two  primary  types  of  undulator  polarization,  linear  and  helical.  In  a  linear 
undulator,  which  is  most  common  in  today’s  FELs,  there  are  two  rows  of  magnets, 
between  which  the  electrons  pass.  These  magnets  alternate  polarities  along  the  undulator, 
so  that  as  the  electron  bunches  pass  along  the  undulator  axis,  they  experience  alternating 
forces,  causing  them  to  weave  with  a  slalom  type  motion.  This  type  of  undulator  is 
shown  in  Figure  2  with  the  electron  path  shown  as  the  light  blue  line. 


Figure  2.  Undulator  and  Resonator  Cavity 

This  slalom  motion  is  caused  by  accelerating  the  electrons  periodically  in  the 
direction  transverse  to  the  motion  of  the  electrons.  In  general,  the  acceleration  of  moving 
charges  in  a  magnetic  field  is  caused  by  the  Lorentz  Force. 
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This  acceleration  causes  the  electrons  to  radiate  along  the  axis  of  motion.  This 
radiation  is  refleeted  off  of  the  mirrors,  whieh  ean  be  seen  in  Figure  2.  The  region 
between  these  mirrors  is  called  the  resonator  eavity.  One  mirror  is  fully  refleeting,  while 
the  other  is  partially  transmissive.  In  the  figure,  the  latter  is  the  mirror  on  the  right.  The 
light  that  is  transmitted  through  the  right  mirror  is  propagated  to  the  target. 

A  helical  undulator  generates  the  magnetie  field  by  wires  wound  around  the 
eentral  eavity  in  a  spiral.  In  this  eonfiguration,  the  eleetrons  see  a  field  whieh  rotates 
around  the  axis  as  they  travel  through  the  undulator.  The  theory  of  this  type  of  undulator 
will  be  examined  in  the  next  seetions. 

The  eleetrons  are  then  deeelerated  through  the  same  aeeelerator  cavities  in  whieh 
they  were  aeeelerated  by  sending  them  through  180°  out  of  phase.  This  offers  two 
advantages.  First,  this  deeeleration  signifieantly  lowers  the  eleetrons’  energy  so  that 
when  they  reach  the  beam  dump,  the  radiation  emitted  by  the  eleetrons  colliding  with  the 
dump  is  significantly  decreased.  Seeond,  this  allows  for  energy  recovery.  The 
decelerating  electrons  transfer  the  energy  that  they  lose  back  to  the  electromagnetic 
fields,  thus  decreasing  the  amount  of  supplied  energy  that  the  fields  need  to  accelerate  the 
electrons.  This  can  significantly  increase  the  wall-plug  efficieney  of  the  FEL,  potentially 
inereasing  a  100  kW  laser  from  2%  to  over  6%.  [12] 

B,  UNDULATOR  FIELDS  AND  THE  RESONANCE  CONDITION 

Within  the  resonator  eavity,  there  are  three  different  eleetromagnetic  fields  whieh 
interaet  with  the  eleetrons.  They  are  the  optieal  magnetic  field,  symbolized  by  Bs,  the 
optical  electric  field,  Es,  and  the  undulator  magnetic  field,  Bu.  These  fields  create  forces 
on  the  electrons  via  the  relativistie  Lorentz  equation,  whieh  in  cgs  units  is 

^(rP)  =  -— (e.+Pxb),  (2.1) 

at  me 

where  v  =  Pc  is  the  relativistic  velocity  of  the  electron,  |e|  is  the  electron  charge 
magnitude,  m  the  eleetron  mass,  c  the  speed  of  light,  B  =  Bs  +  Bu,  and 
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is  defined  as  the  Lorentz  Factor.  The  electron’s  energy  is  then  Ee  =  ymc^,  which  typically 
is  approximately  100  MeV,  yielding  a^^  =  0.999986  «  1,  and  y  w  200. 

In  a  helical  undulator,  the  above  mentioned  fields  are  given  as 


j  =E'(sin^,cos^4/',0) , 

(2.3) 

E'(cos^,-sin^,0) ,  and 

(2.4) 

5  (cosCkgZ),  sin(A:oZ),  O) , 

(2.5) 

where  E  is  the  magnitude  of  the  optical  field,  \j/  =  kz-(X)t  +  (j)  ,k  =  In/X  is  the  optical  wave 
number,  X  is  the  optical  wavelength,  z  is  the  longitudinal  distance  along  the  undulator, 
(u  =  A:c  is  the  optical  frequency,  (j)  is  the  optical  phase,  B  is  the  magnitude  of  the 
undulator  magnetic  field,  and  ko  =  2k/Xq  is  the  undulator  wave  number  with  Xq  the 
undulator  wavelength. 

Resonance  is  an  important  concept  in  the  operation  of  the  free  electron  laser.  It  is 
defined  as  the  conditions  under  which  exactly  one  wavelength  of  light  passes  over  an 
electron  in  one  undulator  period.  See  Figure  3  below,  in  which  the  red  circle  is  as 
electron,  the  blue  line  a  wavelength  of  light,  and  the  green  line  an  undulator  wavelength. 


The  resonance  condition  allows  for  determining  the  wavelength  of  light  that  is 
generated  by  a  FEL.  The  time  that  it  takes  one  electron  to  travel  one  undulator  period  is 
defined  as  At,  and  the  velocity  in  the  z  direction  of  the  electron  is  Vz  =  XqIAI.  In  this  same 
amount  of  time,  the  light  wave  travels  Xq+X  at  speed  c,  which  means  that 

=  ^  =  (2.6) 

Also,  it  has  been  stated  that  v  =  Pc,  and 


9 


p  =  (A, 

from  which 


=c4f-pl  , 


(2.7) 

(2.8) 


(2.9) 


(2.10) 


where  is  the  transverse  veloeity  of  the  eleetron.  It  will  be  shown  in  the 

next  seetion  that  the  eleetrons’  transverse  veloeity  in  the  undulator  is 


P^=  — 

Pl  2  ’ 

Y 


(2.11) 


where  the  dimensionless  undulator  parameter  is  defined  as  if  =  eBXo/2imic^ ,  whieh  is  on 
the  order  of  unity. 

Therefore,  equation  (2.6)  ean  be  written  as: 


Ag  A  -h  A'q 


2  -2 

and  from  equation  (2.2),  P  =  7  -  y'  ,  so 


An 


■  —  A  A, 


'O’ 


An 


■  —  A  -h  A, 


.2 


O’ 


1-r 

r 


and 


An 


■  —  A  +  An 


1- 


1  +  K^ 


(2.12) 


(2.13) 
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Assuming  y»  \  and  K  x\,  then 


1- 


\  +  K^ 


1- 


1  +  K^ 


(2.14) 


Substituting  the  approximation  from  equation  (2.14)  into  (2.13)  yields 

.  .  1  +  ^" 

2  •  (2-15) 

2y 

Equation  (2.15)  yields  the  wavelength  of  the  light  generated  from  the  undulator  at 
resonanee.  It  is  important  to  note  that  the  wavelength  depends  upon  the  strength  of  the 
undulator  magnetie  field  (whieh  is  ineorporated  in  the  dimensionless  undulator  parameter 
K),  the  undulator  period,  and  the  eleetron  energy. 

C.  ELECTRON  MOTION 

Mieroseopie  eleetron  motion  within  the  undulator  and  laser  beam  is  governed  by 
the  pendulum  equation.  Another  eomponent  of  the  relativistie  Lorentz  foree  equation  is 
termed  the  energy  equation,  whieh  is  given  as 

^  =  -  — P-E,  (2.16) 

at  me 


and  is  an  expression  for  the  evolution  of  the  eleetron  energy  over  time  within  the 
undulator.  [9] 

Initially,  it  is  assumed  that  there  is  no  light  in  the  eavity,  so  the  relativistie  Lorentz 
foree  equation  is 


Now, 


(PxB). 

at  me 


pxB  =  5 


X  y 

A  Py 

cos(A:oz)  sin(A:oz) 


(2.17) 


(2.18) 


so  from  equation  (2.17)  and  (2.18)  we  have 
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d_ 

dt 


(A)  = 


eBp, 


-sin(A:oz) 

ymc 
eB 


ymc^kQ 

cBAq 

Inymc^ 


-cos 


[k^z]. 


-cos 


[k^z). 


A  =-— cos(^oz),  (2.19) 

where  K=eBXo/27rmc^  again  is  the  dimensionless  undulator  parameter,  and  the  eonstant  of 
integration  is  zero  for  perfeet  injeetion. 

By  a  similar  argument,  it  ean  be  shown  that 

p  =-—sm{k,z).  (2.20) 

y 

Combining  these  two  equations  yields 

= - (eosA:oZ,sinA:oZ,0) .  (2.21) 

r 

This  gives  the  veloeity  of  the  eleetrons  perpendieular  to  the  longitudinal  direetion. 
The  magnitude  of  P^  is  approximately  0.01,  whieh  shows  that  the  perpendieular  veloeity 
of  the  eleetrons  is  mueh  smaller  than  their  veloeity  along  the  axis. 


Assume  that  the  longitudinal  position  of  the  eleetrons  is  given  by  z(t)  =  f^^ct  ~  ct, 
so  that  A:oz  «  kffit  =  coot.  This  means 

P^  = - (eos(»ot,sin(2;ot,0) .  (2.22) 

r 

Integrating  to  get  the  eleetrons’  transverse  position  yields 

K  X 

= - ^(-sin(2;ot,eosryot,0) ,  (2.23) 

yin: 
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where  eonstants  of  integration  are  taken  to  be  zero  for  perfeet  orbits.  For  a  high  power 
FEL,  y  «100,  and  /lo~2.5  em,  so  that  Xx~40pm.  The  beam  radius  is  typieally 

«40pm. 

Now  light  will  be  ineorporated  to  demonstrate  its  effeet  on  the  mieroseopie 
eleetron  motion.  The  fields  were  given  in  equations  (2.3)  and  (2.4).  Inserting  these  into 
equation  (2.1)  yields 

—  (Py)  =  — —\^E{\-p^^{co?,y/,-?,my/,Q))  +  y0^5(-sin A:QZ,eos A:oZ,0)]  .  (2.24) 

dt  tnc  ' - ^ _ '  's _ ^ _ ^ 

optieal  field  foree  undulator  field  foree 

Using  equations  (2.4),  (2.16),  and  (2.21)  we  find 

/  =  ^^^(cos(kQz)cosi//-sm(kQz)smi//) , 

/me 

pKF 

and  /  = - (eos(A:QZ  +  ^))  (2.25) 

yme 

where  y/  =  kz  -  cot  +  (j) .  Defining  C  =  (k+ko)z-cot  as  the  eleetron  phase  shows  that  equation 
(2.25)  ean  be  rewritten  as 

y  =  ^^cos(<;  +  ^).  (2.26) 

/me 

In  order  for  there  to  be  an  inerease  in  the  optieal  energy,  the  eleetrons  must  lose 
energy.  Therefore,  /  must  be  less  than  zero  when  averaged  over  all  eleetrons  in  the 
beam. 

In  order  to  explore  eleetron  bunehing,  the  eleetron  phase  ^  will  be  examined. 

=  [k  +  kf^)z-cot, 

Q  =  (A:  +  A:q)z-(», 
so  that  Q  -{k  +  ko)  fd/e - oo. 

From  equation  (2.2),  it  ean  be  shown  that 
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Using  the  approximation  from  equation  (2.14), 


(2.27) 


Inserting  this  into  the  above  equation  for  ^  yields 


C  =  {k  +  ko)  1 — c-o). 
V  w  J 


(2.28) 


Since  2<<2o,  then  k»ko,  so  that  ko  is  negligible  compared  to  k.  Making  these 
approximations  and  differentiating  equation  (2.28)  yields 


A  Ir  +  ^  1  • 

^  =  kc  — ^  y. 

y  r  ) 


(2.29) 


Substituting  y  from  equation  (2.26)  into  the  right-hand-side  yields 


^  =  kc 


\  +  K^  1  eKE 


y  j  y  me 


Kf+t*)- 


(2.30) 


From  equation  (2.15)  and  the  fact  that  kc=co=2nc/X,  we  find  that 


E(u¥)’ 


2r" 

kc  =  knC  - ^ 


(2.31) 


Substituting  this  into  equation  (2.30)  yields 


ILeEK 


(2.32) 
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The  concept  of  dimensionless  time  is  introduced  to  further  clarify  equation  (2.32). 
Let  X  =  ct/L,  where  L  is  the  overall  length  of  the  undulator,  so  that  r  =  0  corresponds  to 


the  entrance  of  the  undulator  and  r  =  1  the  exit.  Therefore, 


dr  =  —dt, 

L 

dr^^  dt^^[c  r 


leEKIdk^ 

1  2 
Y  me 


+^^)  • 


(2.33) 


Here  (...)  symbolizes  ^— (...),  the  second  dimensionless  time  derivative.  Now 


L  =  N  A,  =■ 


or  Lkr,  =  ItxN. 


Inserting  this  into  equation  (2.33)  yields 


AttNcKLE 

2  2  ' 

Y  me 


(2.34) 


Defining  the  optical  field  amplitude  as 


ATxNeKL 


«  — —E 


y  me 


(2.35) 


allows  equation  (2.34)  to  be  written  as 


=  a  cos(i^  +  (/)) . 


(2.36) 


This  demonstrates  that  the  electron’s  phase  behavior  within  the  undulator  is 
governed  by  the  pendulum  equation,  driven  by  the  optical  field  of  amplitude  |a| .  The 
dimensionless  phase  velocity  of  the  electron  is  then  defined  as 


■  =  ^  =  L[{k  +  k,)l3^-k] 


(2.37) 


15 


Equations  (2.36)  and  (2.37)  are  significant  because  they  facilitate  the  modeling  of  the 
microscopic  electron  motion  within  the  undulator,  which  will  be  explored  in  a  later 
section. 

D.  THE  OPTICAL  WAVE  EQUATION 

The  final  step  in  describing  the  theory  behind  the  Free  Electron  Easer  is  to 
analyze  the  interaction  between  the  electrons  and  the  optical  wave  they  create.  The 
optical  wave  equation  is 


( 

V 

V 


2 


A(x,0 


4^  X  /  N 

- J^(x,0, 

c 


(2.38) 


in  which  the  vector  potential,  A,  for  a  helical  undulator  is  given  by 


(sin^,cos^i^,0) 


where  y/  =  kz  -  cot (f) ,  from  which 


(2.39) 


=VxA 


(2.40) 


and 


c  dt 


(2.41) 


yields  the  specified  field  Bs  and  Es,  and  is  the  transverse  current  density. 


In  order  to  solve  this  equation,  a  number  of  assumptions  will  be  made.  First,  the 
optical  waves  are  assumed  to  be  plane  waves,  so  that  there  is  no  change  in  the 
perpendicular  components  of  A,  and  only  the  longitudinal  components  of  the  Eaplacian 
survive.  This  means  that 


V"A 


(2.42) 


The  second  assumption  is  that  the  fields  are  slowly  varying  in  z  and  t.  This  means  that,  in 
the  time  domain,  during  an  optical  period  E  <sc  coE  and  ^  <sc  ax!) ,  and  in  space,  over  an 
optical  wavelength  E"  and  .  [3  p.l3]  This  allows  all  second-order  and 
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higher  terms  to  be  neglected.  When  equation  (2.39)  is  substituted  into  (2.38),  these 
assumptions  simplify  the  expression  to 


^ BE  ^  \  BE' 
\Bz  cBt 


Bcj)  1  Bcj)" 


An 


(cos^i^,-sin^,0)-2£'  — H - (sin^,cos^i^,0)  « - Jj^  .(2.43) 

\Bz  c  Bt )  c 


In  order  to  express  equation  (2.43)  in  dimensionless  parameters,  the  spatial  and 
time  derivatives  can  be  expressed  in  terms  of  the  undulator  length  L  and  dimensionless 
time  T.  The  spatial  coordinate  z*  is  used  to  define  a  point  that  follows  the  light  and  is 
defined  as 


z*=z-ct.  (2.44) 

From  the  definition  of  r  =  ct/L,  equation  (2.44)  can  be  rewritten  as 


z*  =z-Lt 

(2.45) 

This  allows  the  partial  derivatives  to  be  written  as 

B  Bz*  B  Bt  B 

Bz  Bz  Bz*  Bz  Bt 

(2.46) 

and 

B  Bz*  B  Bt  B 

Bt  Bt  Bz*  Bt  Bt 

(2.47) 

Using  the  definition  of  r  and  equation  (2.45),  5z’  /  5z  =  1 ,  Bz*  /  Bt  = 
Bt / Bt  =  c/ L  .  Equations  (2.46)  and  (2.47)  then  become 

-c ,  Bt ! Bz  =  0,  and 

B  _  B 

Bz  Bz* 

(2.48) 

and 

B  B  c  B 

—  —  ~C  — H - . 

Bt  Bz  L  Bt 

(2.49) 

Combining  equations  (2.48)  and  (2.49)  yields  the  operator 


5^15  ^  ^ 

Bz  c  Bt  Bz*  c 


B  c  B 

~C - H - 

Bz  L  Bt 


lA 

L  Bt 


(2.50) 


When  this  operator  is  applied  to  equation  (2.43),  it  is  then  written  as 


17 


2  BE ,  .  .  2E  d(p  i  .  An 

L  dr  L  dr  c 


(2.51) 


In  order  to  solve  this  equation,  the  veetor  eomponents  must  be  deeoupled.  This  is 
done  by  defining  two  unit  veetors,  =  (oos^,-sin^^,0)  and ^2  =  • 

Multiplying  equation  (2.51)  by  these  two  unit  vectors  individually  yields 


BE  2nL  , 

—  = - Ji 

dt  c  ^  ' 


(2.52) 


and 


j^d<f)  _  2nL 
dt  c 


J±*^2  • 


(2.53) 


The  transverse  current  density  is  generated  by  the  transverse  motion  of  the 
electron  within  the  undulator,  and  is  the  sum  of  the  individual  electron  currents,  given  by 

J,=2;[-ecP,5>(x-r,)],  (2.54) 

i 

where  (5^(x-ri)  is  the  three  dimensional  Dirac  delta  function  and  ri  is  the  position  of  the  i*’’ 
electron.  Inserting  this  expression,  as  well  as  Pj^  from  equation  (2.21),  yields 


dt 


2nKeL 

r 


^  (5'^  (x  -  r^. )  (cos(A:oz)  cos  y/  -  sin(A:Qz)  sin  y/^ 


and 


(2.55) 


E  —  =  V  - r. ) (cos(A:oz) sin  yr  -  sin^k^z)  cos .  (2.56) 

dr  y 

The  above  summation  of  electron  phases  can  be  replaced  by  the  volume  density 
of  electrons  and  an  average  over  sampled  electrons  <...>.  This,  along  with  trigonometric 
identities  and  the  definition  4”=  {k+kk)z-cot  given  in  section  2.C  yields 

8E_2n^  (2.57) 

5r  y 

and 
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(2.58) 


^8l^2n^ 

dry 

where  <...>  is  the  average  over  all  sampled  eleetrons.  Multiplying  equation  (2.58)  by  i 
and  adding  it  to  (2.57)  reeombines  the  separated  equations,  and  yields 


dr 


dr 


r 


(2.59) 


Multiplying  both  sides  of  the  equation  by  allows  it  to  be  written  as 


dr 


InKeLp^ 

r 


(2.60) 


Equation  (2.35)  ean  be  written  as 


2  2 

a  y  me 

AnNeKEL 


(2.61) 


so  the  left  side  of  equation  (2.60)  ean  be  multiplied  by  (2.61)  without  loss  of  equality,  and 
ean  therefore  be  written  as 


(2.62) 


where  a  =  |a|e'^  is  the  eomplex  dimensionless  optieal  field  and  j  is  defined  as  the 
dimensionless  eurrent  density,  whieh  is 


j  = 


3  2 

Y  me 


(2.63) 


Equation  (2.62)  demonstrates  that  the  evolution  of  the  optieal  field  is  dependent  upon  the 
eleetron  beam  eurrent  and  the  average  phase  of  the  eleetrons. 

E,  GAIN 

Gain,  G,  is  defined  as  the  fraetional  ehange  of  power  in  the  optieal  field  per  pass 
through  the  undulator,  and  is  given  by 

G  =  (2.64) 
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where  ao  and  a\  are  the  dimensionless  optieal  field  strengths  at  r  =  0  and  r  =  1,  i.e.  the 
beginning  and  end  of  the  undulator,  respeetively.  In  eaeh  pass  through  the  undulator, 
energy  is  transferred  between  the  eleetrons  and  the  optieal  field.  Therefore,  the  gain  ean 
be  analyzed  by  determining  the  ehange  in  the  energy  of  the  eleetrons. 

It  will  first  be  shown  that  a  ehange  in  the  eleetron  energy  is  proportional  to  its 
ehange  in  phase  veloeity.  This  phase  veloeity  v  was  given  in  equation  (2.37)  as 

v  =  L\{k  +  k^)(i^-k].  (2.65) 

Beeause  ko«k  for  relativistie  eleetrons,  this  ean  be  written  as 

A.V  X  LkA./3^  (2.66) 


Using  k  =  2n  / X  and  the  resonanee  eondition  defined  in  equation  (2.15),  with  N  =  L  / Xq, 
this  ean  be  written  as 


Av  » InN 


AA 


In  seetion  2.C,  it  was  shown  that 


Pl=^- 


\  +  K^ 


Taking  the  time  derivative  of  both  sides  and  solving  for  yields 


y  \  +  K^ 

7~^ 


(2.67) 


(2.68) 


(2.69) 


Multiplying  both  sides  by  At  and  using  lS./3^  =  y^^At  and  A/^  =  /^At  yields 

Ay(l  +  .^") 


Ay0.=- 


(2.70) 


Using  the  approximation  that  Pz~\,  equation  (2.70)  is  inserted  into  equation 
(2.67)  to  yield  the  relationship  between  Av  and  Ay,  whieh  is 


Av  =  4;riV^. 

r 


(2.71) 
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Solving  for  Ay  and  multiplying  both  sides  of  the  resulting  equation  by  me  yields 


(Ay)mc^ 


yme^ 

4/rN 


(2.72) 


where  <v>  is  the  average  eleetron  phase  veloeity  for  all  eleetrons  at  the  end  of  the 
undulator  and  vo  is  the  initial  eleetron  phase  veloeity.  Sinee  ymc  is  the  energy  of  the 
eleetrons,  this  equation  shows  the  ehange  in  an  individual  eleetron’s  energy  in  terms  of 
the  ehange  in  eleetron  phase  veloeity.  The  number  of  eleetrons  within  a  small  a  length 
dV  in  an  optieal  wave  is  dNe  =  pFdV  where  F  is  the  fdling  faetor,  defined  as  the  eross- 
seetional  area  of  the  eleetron  beam  divided  by  the  eross-seetional  area  of  the  optieal 
beam.  The  initial  energy  of  the  volume  element  of  the  optieal  beam  is 


E 

dEr.  = — dV , 
°  4;r 


(2.73) 


for  eireularly  polarized  light.  The  gain  is  ealeulated  assuming  the  lost  energy  of  the 
eleetron  beam  is  eonverted  to  laser  energy,  so  that 


AE 

^  _  electron-beam 


pEdV 


ymc^[iy)-v^) 


4nN 


dE, 


optical-beam 


An 


dV 


pEymc\{v)-v,) 

Lr  = - z - 

NE^ 

From  equations  (2.35)  and  (2.63) 

j  _  ypme^ 
al  ~  2NE^  ■ 


(2.74) 


(2.75) 


(2.76) 


Substituting  this  into  equation  (2.75)  allows  the  gain  equation  to  be  written  as 

G  =  -^-({v)-v,).  (2.77) 

a. 

This  equation  shows  that  in  order  for  there  to  be  positive  gain,  the  average  phase  veloeity 
must  be  less  than  the  initial  eleetron  phase  veloeity.  Therefore,  there  has  to  be  a  net  loss 
of  eleetron  phase  veloeity  in  order  to  aehieve  positive  gain  in  a  FEL. 
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F.  PHASE  SPACE 

Phase  space  is  a  useful  graphical  technique  for  understanding  the  physics  of  the 
pendulum  equation.  As  was  shown  in  section  2.C,  electron  motion  within  the  undulator 
is  governed  by  the  pendulum  equation  (equation  (2.36)),  so  the  use  of  phase  space  is  an 
ideal  method  to  show  the  evolution  of  the  electrons.  The  electrons’  phase  ^  and  phase 

velocity  v  =  (^  are  plotted  on  the  horizontal  and  vertical  axes,  respectively,  as  in  Figure  4 
below,  over  one  pass  through  the  undulator. 

B 


-5 

-^1/2  ^  3n/2  0  T  L 

Figure  4.  Phase  Space  Plot  for  Electrons  on  Resonance 
The  width  of  the  window  is  a  section  of  the  electron  beam  one  wavelength  of  light 
long,  and  moves  along  at  the  speed  of  an  electron  at  resonance.  In  Figure  4,  20  sample 
electrons  are  followed  from  r  =  0  to  r  =  1 ,  symbolized  by  the  change  of  the  color  of  the 
tracks  from  yellow  to  red,  respectively.  Because  the  electrons’  energy  change  is 
proportional  to  v,  this  plot  demonstrates  the  evolution  of  the  electrons’  energy  as  well.  In 
this  case,  all  electrons  start  at  an  initial  phase  velocity  of  vo  =  0,  or  at  resonance,  and  are 
uniformly  distributed  from  -%I2  <  C  <  37i/2.  It  can  be  seen  in  Figure  4  that  while  the 
electrons  that  begin  below  C  =  0  gain  energy,  those  greater  than  C  -  0  lose  the  same 
amount.  This  means  that  there  is  no  net  transfer  of  energy  from  the  electrons  to  the 
optical  field. 
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The  red  line  around  the  eleetrons  is  the  seperatrix,  which  shows  the  demarcation 
between  open  and  closed  orbits.  The  electrons  that  start  outside  this  line  will  have  open 
orbits,  while  those  that  are  contained  within  it  follow  closed  orbits. 

In  Figure  5  below,  the  electrons  are  injected  with  an  initial  phase  velocity  vq  =  2.6, 
again  uniformly  distributed  along  f. 


Figure  5.  Phase  Spaee  Plot  With  Vo=2.6. 

Note  that  in  Figure  5,  more  electrons  lose  energy  than  gain  energy.  This  net  loss 
of  electron  energy  means  that  there  is  gain  in  the  optical  field.  This  gain  is  graphed  in  the 
upper  right  hand  side  of  the  figure. 

Figures  4  and  5  are  idealized  examples  of  simulation  outputs.  Real  electrons  will 
not  all  arrive  at  the  undulator  with  the  same  Vo,  so  electrons  within  simulations  are  given 
an  initial  energy  spread  in  order  to  accommodate  this  practical  issue.  The  use  of  sample 
electrons  increases  the  speed  of  simulations,  but  may  decrease  accuracy.  Therefore  many 
more  sample  electrons  than  the  20  shown  in  the  above  figures  are  used.  Simulations 
conducted  in  order  to  optimize  FEL  performance  use  this  and  other  techniques  to 
demonstrate  the  output  that  can  be  achieved  for  a  given  set  of  parameters.  This  will  be 
further  explored  in  later  sections. 
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III.  MULTIPLE  BEAM  DIRECTORS 


The  theory  whieh  deseribes  a  megawatt-elass  free  eleetron  laser  is  essentially  the 
same  as  that  for  existing  FELs.  Ongoing  engineering  improvements  opties  and  eathodes 
indieate  that  megawatt  operation  ean  be  realized  in  the  future. 

However,  one  thing  that  will  limit  all  weapons-grade  lasers  is  thermal  blooming. 
Thermal  blooming  is  eaused  by  heating  of  the  atmosphere  as  a  result  of  absorption  of 
laser  beam  energy  traveling  through  a  volume  of  air.  Beeause  the  optieal  beam  has  a 
Gaussian  eross  seetion,  with  its  peak  in  the  eenter,  the  heating  is  more  signifieant  in  the 
eenter  of  the  beam  than  the  outer  edges.  The  resulting  gradient  in  the  refraetive  index  of 
the  air  through  whieh  it  is  passing  eauses  the  light  to  bend  away  from  the  warmer  eenter 
of  the  beam,  and  thus  behave  as  if  it  were  going  through  a  diverging  lens.  This 
defoeusing  ean  greatly  degrade  laser  effeetiveness. 

As  has  already  been  mentioned,  the  free  eleetron  laser  possesses  the  important 
attribute  of  tunability.  This  is  important  beeause  there  have  been  numerous  studies  that 
investigate  whieh  wavelengths  of  light  pass  through  the  atmosphere  with  the  least 
absorption.  Figure  6  below  shows  the  pereent  transmission  for  light  of  various 
wavelengths.  It  also  shows  at  whieh  wavelengths  absorption  by  CO2  and  H2O  moleeules 
is  most  signifieant. 


Figure  6.  Transmittance  of  1000  ft  horizontal  path  at  sea  level  at  79°F  containing  5.7nm  precipitable  water 

tl3] 
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Because  of  tunability,  the  PEL  can  produce  light  at  the  best  wavelength  for 
maritime  applications,  which  has  been  determined  to  be  approximately  1  pm.  Despite 
this  advantage,  thermal  blooming  is  still  a  significant  problem.  One  way  to  combat  this  is 
the  use  of  multiple  beam  directors,  which  is  the  topic  of  this  chapter. 

A.  DISCUSSION 

There  are  numerous  advantages  to  using  multiple  beam  directors.  By  sending  a 
fraction  of  the  total  energy  through  each  volume  of  space,  the  heating  of  the  atmosphere 
is  reduced.  Therefore,  thermal  blooming  can  occur  far  enough  along  in  the  beam  path  as 
to  have  little  to  no  impact,  which  is  the  topic  of  the  next  section.  This  can  allow  for 
greater  power  on  target  because  these  beams  can  be  aimed  such  that  they  are  focused  on 
the  same  spot.  While  less  energy  passes  through  each  volume  of  air,  the  total  energy 
desired  is  still  focused  to  one  spot  at  the  target. 

Rather  than  combining  the  beams  at  the  target,  another  possibility  is  to  focus  each 
director  on  a  different  section  of  the  target.  For  targets  of  sufficient  size,  directors  could 
be  focused  on  different  sections  in  order  to  gain  the  advantages  of  each  aim  point.  For 
example,  one  could  be  focused  on  the  nose  cone  in  order  to  destroy  the  aerodynamics, 
guidance  systems,  and  other  components  contained  within,  while  the  other  focuses  on  a 
flight  control  surface  such  as  a  wing  to  destroy  the  stability  and  control.  This  would 
increase  the  dwell  time  required,  but  may  be  a  viable  option  for  lower  velocity  targets. 
Another  option  for  targets  that  allow  a  longer  dwell  time  is  the  ability  to  make  the 
affected  area  bigger.  Instead  of  focusing  the  beams  on  the  same  spot,  they  can  be  focused 
adjacent  to  each  other,  doubling  the  affected  area. 

One  drawback  of  this  may  be  the  additional  space  required  on  the  ship.  Directors 
are  large  and  complicated,  as  can  be  seen  in  a  conceptual  drawing.  Figure  7  below. 
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Figure  7.  Conceptual  laser  beam  director. 

This  complication  has  already  been  addressed,  however.  In  most  proposals  for 
laser  placement  on  ships,  two  directors  are  required  in  order  to  ensure  360°  coverage. 
Figure  8  below  shows  one  such  proposal  for  fore-and-aft  mounted  directors.  The  left 
picture  shows  a  close-up  view  of  the  firing  arcs,  while  the  right  is  expanded  to  show  the 
overlap  areas  where  two  beams  can  hit  the  same  target.  This  shows  that  the  angular 
coverage  of  which  two  directors  are  capable  is  approximately  80°  on  each  side,  allowing 
for  the  beams  to  be  able  combine  on  a  single  target  for  almost  half  of  the  total  coverage 
area.  This  is  by  no  means  the  only  configuration  that  can  be  used  in  order  to  achieve 
360°  coverage,  but  it  is  the  most  common.  Because  two  beam  directors  will  most  likely 
be  placed  on  the  ship  for  this  reason,  it  makes  sense  to  develop  the  necessary  techniques 
in  order  to  direct  both  beams  to  a  single  target  in  any  direction.  The  advantages  of 
different  placement  options  will  be  explored  in  a  later  section. 
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B, 


EFFECT  OF  DISTANCE  ALONG  BEAM  OF  THERMAL  BLOOMING 


1,  Simulation  Methods 

As  was  mentioned  above,  the  most  signifieant  way  that  using  multiple  beam 
direetors  will  help  deerease  thermal  blooming  is  the  faet  that  the  point  where  it  takes 
plaee,  ealled  the  stagnation  zone,  is  moved  to  near  the  end  of  the  beam’s  path.  This 
allows  less  time  for  the  beam  to  diffraet,  so  the  spot  size  is  not  inereased  as  mueh  as  it 
would  be  if  thermal  blooming  took  effeet  earlier  in  the  beam. 

In  order  to  study  this,  eomputer  models  of  the  beam  path  were  used  to 
demonstrate  the  effeet  of  moving  the  stagnation  zone.  This  simulation  method  is  shown 
below  from  [5]. 


The  parabolie  wave  equation  deseribes  the  propagation  of  a  laser  beam  to  a  target; 


5 


(3.1) 


where  the  eomplex  field  a  =  |a|e'^  evolves  in  dimensionless  time  Xr  =  0^1  along  the 

range  to  target,  and  =  5^  +  5^  is  the  dimensionless  transverse  Laplaeian.  Equation 

(3.1)  uses  the  dimensionless  field  amplitude |a| ,  whieh  is  proportional  to  the  aetual  optieal 

field  amplitude  [eI  ,  as  derived  in  Chapter  II.  The  initial  Gaussian  field  shape  at  =  0  is 

given  hy  a{r)  =  aQe  ''  where  =  ,  r'  =  x  ^  y^,  ao  =  1,  and 

=  1  +  /  Zq  ,  for  a  dimensionless  Rayleigh  length  zo  and  beam  waist  at  r^- 

The  result  of  the  simulation  yields  a  matrix  of  optieal  amplitudes,  whieh  are  then 
plotted.  Figure  9  below  shows  the  side  and  end-on  views  of  a  beam  traveling  through  the 
atmosphere  with  no  blooming  effeets. 
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Figure  9.  Beam  propagation  with  no  thermal  blooming. 


The  color  scale  on  the  right  of  the  figure  shows  the  amplitude  of  the  optical  field, 
blue  being  the  lowest  and  red  the  highest.  The  left  side  of  the  figure  shows  a  side  view  of 
the  beam  as  it  travels  from  a  source  on  the  left  (x^  =  0)  to  target  on  the  right  (x;-  =  1).  It  is 
expanded  for  the  purposes  of  illustration,  as  the  actual  beam  diameter  tapers  from  a  meter 
to  five  centimeters,  and  the  distance  is  on  the  order  of  eight  kilometers.  The  right  side  is 
an  end-on  view  of  the  beam  at  full  range.  Notice  the  clearly  defined,  concentric  circles  of 
intensity.  This  is  due  to  the  Gaussian  intensity  profile  mentioned  in  the  introduction  of 
this  section.  The  area  of  this  circle  is  the  focus  of  this  analysis. 

The  stagnation  zone  is  modeled  by  inserting  a  phase  distortion  into  the  path  of  the 

beam. 


(3.2) 


where  dcj)^  is  the  effective  zone  strength  of  the  lens. 

2.  Simulation  Results 

The  next  figure  shows,  from  top  to  bottom  the  effect  of  placing  a  stagnation  zone 
at  distances  of  30%,  50%,  and  80%  (x,-  =  0.3,  0.5,  and  0.80,  respectively)  along  the  beam 


path. 
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Figure  10.  Stagnation  zone  at  (a)  30%,  (b)  50%  and  (c)  80%  of  total  range. 

Notice  the  diffraction  of  the  beam  and  the  inereased  size  of  the  spot  on  the  right  side  of 
the  figures.  The  total  area  deereases  as  the  zone  is  moved  further  along  the  beam. 

The  strength  of  the  stagnation  zone  in  the  model  ean  be  varied  as  well  as  the 
loeation.  In  Figure  1 1  below,  the  area  of  the  spot  is  plotted  as  a  funetion  of  the  distanee 
along  the  beam  for  various  zone  strengths.  The  vertical  axis  is  normalized  to  the  radius 
of  a  spot  that  is  unaffeeted  by  thermal  blooming,  where  A  is  the  area  of  the  affeeted  spot, 
and  Ao  the  area  of  the  unaffeeted  spot.  The  horizontal  axis  is  Tc,  the  loeation  along  the 
dimensionless  path  where  the  stagnation  zone  is  loeated.  The  strength  of  the  zone  is 
noted  in  the  legend,  which  is  the  value  of  from  equation  (3.2). 
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Figure  11.  Positioning  of  the  stagnation  zone  for  varying  strength  of  the  zone. 

The  area  of  the  spot  for  beams  affeeted  by  blooming  is  much  larger  than  the  area 
of  a  spot  from  an  unaffected  beam.  It  reaches  a  maximum  when  blooming  occurs  at 
approximately  55%  of  the  range  to  target,  and  then  drops  sharply  after.  An  interesting 
feature  of  this  plot  is  the  aforementioned  peak.  This  is  due  to  the  area  of  the  beam  at  the 
time  blooming  takes  effect.  Near  the  source  of  the  beam  at  the  ship,  the  beam  has  a 
relatively  large  area,  and  therefore  the  intensity  is  lower  so  blooming  tends  not  to  occur. 
As  can  be  seen  in  equation  (3.2),  the  strength  of  the  stagnation  zone  is  dependent  on  the 
amplitude  of  the  optical  field,  so  the  lens  effect  is  weaker  than  it  is  later  in  the  beam.  The 
fall  off  of  the  area  following  the  maximum  is  due  to  the  lack  of  distance  after  blooming 
occurs  for  the  beam  to  expand.  Even  though  the  blooming  is  stronger,  as  the  stagnation 
zone  moves  toward  the  target,  the  beam  has  less  time  to  expand,  so  the  spot  size 
decreases  as  the  zone  approaches  the  termination  of  the  beam.  These  effects  can  also  be 
seen  in  Figure  10. 

Figure  12  shows  a  more  complete  three  dimensional  plot  of  the  effect  of  the 
strength  and  position  of  the  stagnation  zone. 
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Figure  12.  3-D  plot  of  stagnation  zone  plaeement  and  strength. 


Both  figures  shows  the  advantage  of  moving  the  stagnation  zone  to  as  elose  to  the 
target  as  possible.  The  spot  size  is  close  to  that  of  an  unaffected  beam  when  the  location 
of  the  stagnation  zone  is  >90%  of  the  distance  to  the  target.  As  will  be  shown  later,  the 
use  of  multiple  beams  causes  the  stagnation  zone,  if  it  does  occur,  to  be  located  at  >99% 
of  the  distance  to  the  target,  so  that  the  area  will  be  barely  affected  by  blooming. 


C.  GEOMETRIC  DERIVATION  OF  BEAM  OVERLAP 


In  order  to  determine  the  advantages  of  multiple  beam  directors,  it  is  necessary  to 
investigate  the  amount  of  overlap  that  the  two  beams  will  experience  in  their  path  to  the 
target.  This  will  allow  for  determining  how  much  the  atmosphere  will  be  heated,  and 
therefore  the  degree  to  which  the  probability  of  thermal  blooming  will  be  lessened. 


This  section  of  the  thesis  will  show  how  the  amount  of  overlap  the  beams  will 
experience  is  determined  as  a  function  of  their  distance  apart  from  each  other  on  the  ship. 
For  the  purposes  of  this  analysis,  the  beams  will  be  assumed  to  be  firing  at  a  spot  that  is 
equidistant  between  the  two  directors.  For  example,  directors  placed  fore-and-aft  will  be 
firing  at  a  target  approaching  directly  perpendicular  to  the  middle  of  the  ship. 
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A  simple  geometric  derivation  will  show  how  this  is  done.  Figure  13  below 
shows  the  termination  of  the  beams,  where  they  overlap  in  which  rt  is  the  radius  of  the 
beam  at  the  target,  /o  is  the  length  of  the  overlap,  /p  is  the  length  before  the  overlap,  A  is 
the  distance  to  the  target,  Xs  is  the  distance  off-axis  of  the  directors,  and  is  the  radius  of 
the  directors  (source). 

n 


Figure  13.  Overlap  of  beams  at  target 

In  order  to  determine  the  percent  of  the  range  that  the  beams  are  overlapping,  look 


at  one  side  of  this  figure. 


rt 
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Figure  14.  Triangular  representation  of  beam  overlap 

Similar  triangles  yield  that 
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(3.3) 


r,  _  k 

Ip 

The  desired  radius  of  the  beam  at  the  target,  r^,  is  about  5  em,  and  the  radius  of  the 
direetor  is  approximately  0.5  m.  Therefore,  given  a  distanee  between  the  direetors,  2xs,  a 
ratio  is  found  using  equation  3.1  that  gives  the  pereentage  of  the  range  over  whieh  the 
beams  are  overlapping,  whieh  will  be  analyzed  in  the  next  seetion.  It  is  worth  noting  that 
this  range  gives  the  initial  point  where  the  beams  start  touehing  eaeh  other.  This  means 
that  the  power  is  not  automatieally  double  at  this  point.  Only  at  the  aetual  aim  point  are 
the  beams  fully  overlapping,  and  doubled  in  power. 

In  order  to  ealeulate  the  volume  that  is  eontained  within  the  overlapping  beams, 
make  the  simplifying  assumption  that  they  are  reetangular  solids.  The  volume  is 


V  =  ^{2r,)l^h  (3.4) 

where  h  is  the  height  of  the  beams,  as  in  Figure  15  below.  The  height  h  is  assumed  to  be 
eonstant,  and  is  estimated  to  be  10  em,  eonsistent  with  the  above  assumption  of  a  5  em 
radius  spot  size. 

h 


Figure  15.  Overlapping  volume. 

D.  PLACEMENT  OPTIONS  FOR  TWO  DIRECTORS 

As  has  been  mentioned,  there  are  numerous  options  for  the  plaeement  of  multiple 
direetors.  The  advantage  of  the  fore-and-aft  plaeement  is  360°  eoverage,  but  this  may 
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also  be  accomplished  by  placing  them  amidships  on  the  port  and  starboard  sides.  Or  they 
could  be  placed  on  the  forward  superstructure  to  one  side,  and  aft  on  the  other.  These 
types  of  configurations  would  allow  360°  coverage,  but  do  not  allow  for  overlap,  so  these 
are  not  ideal  cases  for  this  proposal.  Therefore,  discussion  will  be  focused  on  fore-and- 
aft  placement,  with  some  consideration  given  to  tandem-mounted  directors. 

Figure  16  below  shows  the  overlap  length  of  the  beams  as  a  function  of  the 
horizontal  separation  of  the  beam  directors.  This  is  given  as  the  percentage  of  the  total 
range  to  the  target  that  the  beams  initially  touch. 


Figure  16.  Percent  overlap  of  two  beam  directors. 

This  shows  that  the  percent  of  the  overlap  decreases  quickly  as  the  directors  are 
separated.  Fore-and-aft  placement  would  extend  the  axis  to  much  larger  values.  The 
typical  length  for  a  cruiser  or  destroyer  is  approximately  155  m.  If  the  directors  were  at 
this  separation,  the  overlap  length  would  be  only  0.03%.  The  length  of  an  aircraft  carrier 
is  approximately  335  m,  which  yields  0.015%.  For  a  range  of  8  km,  this  yields  a  distance 
from  initial  contact  of  the  two  beams  to  termination  at  the  target  of  only  2.6  m  and  1.2  m, 
respectively. 

If  it  were  determined  that  placing  two  directors  on  one  mount  was  a  desirable 
alternative  in  order  to  take  advantage  of  multiple  beams,  then  they  could  conceivably  be 
placed  approximately  2  m  apart.  From  the  above  graph,  this  yields  2.5%,  which  means 
an  overlap  distance  of  only  200  m  at  8  km  range. 
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Using  the  geometrical  derivation  in  section  3.B,  an  approximate  volume  of  the 
overlapping  sections  can  be  determined.  This  discussion  will  look  at  the  155m  and  335m 
separation  as  shown  in  the  following  graph. 


range  [m] 


Figure  17.  Overlap  volume  for  two  direetors. 

The  largest  volume  is  obviously  at  the  longest  range,  where  there  is  the  most 
overlap.  However,  even  this  number  is  very  small,  0.029  m  for  155  m  separation,  and 
0.013  m^  for  335  m. 


As  a  comparison,  the  overall  volume  that  one  beam  contains  will  be  calculated. 
The  formula  for  the  volume  of  a  tapered  cylinder  is 


V  = 


nl,[ 


2  2 
+r/^+r, 


(3.5) 


where  k  is  the  range  to  target,  is  the  diameter  of  the  beam  director,  and  rt  is  the  beam 
radius  at  the  target.  Inserting  the  assumed  values  of  k  =  8000m,  =  0.5m,  and  rt=  0.05m 

yields  an  overall  volume  of  2320  m  .  This  means  that  the  overlapped  volume  is  a  very 
small  fraction  of  the  overall  volume  of  the  beams. 
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Another  consideration  is  the  motion  of  the  target.  As  was  stated  above,  for  fore- 
and-aft  directors,  there  are  only  a  few  meters  over  which  the  beams  overlap.  Because  of 
the  motion  of  the  target,  however,  this  overlap  does  not  stay  in  the  same  volume  of  air,  so 
it  will  not  heat  sufficiently  for  blooming  to  occur.  While  the  target  moves,  the  air  in  the 
volume  of  the  combined  beams  is  constantly  changing.  A  modest  velocity  for  a  cruise 
missile  is  approximately  mach  0.8,  or  around  270  m/s.  Obviously,  it  will  cover  this 
overlapped  volume  in  a  small  fraction  of  a  second,  so  there  is  little  time  for  the  air  to  heat 
enough  to  cause  thermal  blooming.  Even  for  beams  originating  2  m  apart,  the  length  of 
the  overlapped  portion  of  the  beam  path  is  on  the  order  of  200  m  at  8  km,  so  the  volume 
of  air  where  the  beams  are  combined  is  refreshed  in  less  than  a  second,  which  is  still  too 
short  a  time  for  thermal  blooming  to  take  effect.  [14] 
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IV.  SIMULATIONS  OF  THE  SHORT  RAYLEIGH  LENGTH  EEL 


Computer  simulations  are  powerful  tools  for  free  eleetron  laser  developers.  They 
allow  researchers  to  vary  parameters  (e.g.  electron  beam  current  and  energy,  undulator 
wavelength,  etc.)  in  order  to  optimize  performance,  understand  operation,  and  determine 
what  changes  can  be  made  in  order  to  increase  extraction  and  power,  prior  to  actual 
construction.  Also,  they  are  helpful  in  showing  the  effects  of  various  external  stimuli  on 
laser  operation,  such  as  mirror  tilt  and  vibration.  In  this  case,  three  dimensional 
simulations  were  used  to  study  the  effects  of  changing  the  number  of  undulator  periods, 
the  electron  beam  focus,  the  normalized  Rayleigh  length,  and  the  mirror  output  coupling, 
as  well  as  the  effects  of  tilted  resonator  cavity  mirrors. 

The  Rayleigh  length  is  the  characteristic  distance  over  which  the  optical  beam 
expands  in  area.  Thus,  shorter  Rayleigh  length  means  that  the  beam  expands  faster. 
Cavity  mirrors  can  reflect  only  a  certain  amount  of  intensity  before  damage  occurs.  By 
using  a  short  Rayleigh  length,  the  mirrors  can  be  closer  together  because  the  beam 
expands  more  quickly,  so  the  power  is  spread  over  a  greater  area  on  the  cavity  mirrors  in 
order  to  decrease  the  intensity.  This  is  important  because  of  the  size  limitations  inherent 
in  placing  a  PEL  on  a  ship.  Also,  the  focusing  allows  for  amplification  of  a  single  optical 
wavefront,  which  greatly  increases  beam  quality.  [4]  In  Figure  18  below,  the  optical 
beam  is  outlined  in  blue,  the  electron  beam  in  red,  the  undulator  in  green,  and  the 
magnetic  field  in  black.  The  electron  beam  can  also  be  focused  in  order  to  optimize  the 
transfer  of  energy  between  the  electrons  and  optical  field. 


L 
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Figure  18.  Optical  and  electron  beam  focusing  in  undulator. 

There  are  two  types  of  simulations  that  will  be  discussed,  single-mode  and 
multimode.  For  both  types,  the  following  parameters  were  used: 
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Undulator  Parameters 

Number  of  Periods 

N=  14 

Length 

E  =  37  cm 

Wavelength 

ko  =  2.7  cm 

Magnetic  Field 

B  =  0.8T 

Undulator  Parameter 

K=  1.4 

Electron  Beam 

Energy 

Ee=100  MeV 

Peak  Current 

Ipeak  =  1500  A 

Average  Current 

lavg  ~  1 . 1  A 

Bunch  Eength 

lb  =  0.3  mm  (1  ps) 

Beam  Eocal  Radius 

rb  =  0.07  mm 

Emittance 

£n  =  YorbOb  =  9mm/mrad 

Eocal  Position 

Zb  =  0.5E 

Angular  Spread 

0x=0y=0z=(0x^+0y^)^'^=O.65mrad 

Optical  Mode 

Wavelength 

k  =  1  pm 

Resonator  Eength 

S  ~  16  m  -  20  m 

Waist  Radius 

Wo  =  0.09  mm 

Spot  Radius  at  Mirrors 

w  =  3  cm 

Table  1.  Simulation  parameters. 


A.  SINGLE-MODE  SIMULATIONS 

The  equations  for  the  averaged  eleetron’s  motion  are  derived  assuming  that 
betatron  focusing  is  small,  and  averaging  over  an  undulator  period  to  remove  the  fast 
wiggling  motion  [6] 


x(z)  =  Xo+^^(z-z,) 

(4.1) 

and 

y(^)  =  yo+^y(^-^b)’ 

(4.2) 

where  0x  and  0y  are  the  injection  angles  for  the  sample  electron,  xo  and  yo  are  the  initial 
offsets,  and  Zb  is  the  location  of  the  electron  beam  waist.  This  is  illustrated  by  the  red 
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lines  in  Figure  18  above,  where  not  all  of  the  eleetrons  are  injeeted  on  axis  or  at  the  same 
angle.  The  variable  z  is  again  the  loeation  along  the  undulator,  as  in  Chapter  II. 

The  optieal  field  is  assumed  to  be  traveling  in  the  z-direetion  with  a  Gaussian 
transverse  profile,  given  by 


E{r,z)  =  E^ 


A 


where 


=  -  tan 


+  - 


7rr-{z-z„) 
Az^ 


is  the  optieal  phase  and 


A  =  /iZg 


1+ 


(4.3) 

(4.4) 


(4.5) 


is  the  optieal  mode  area. 

The  evolution  of  the  eleetron  energy  Ee  =  ymc^  along  the  undulator  is  given  by 


'eK\E,\' 

^  ymc^  j 


[Jo(^)-J,(^)]eos(^  +  ^) 


(4.6) 


2  2 

where  Jo  and  Ji  are  Bessel  funetions,  ^  =  K  /2(1  +K  ),  and  C  =  (k+ko)z-a)t  is  the  eleetron 
phase  as  in  Chapter  II. 


Integrating  the  eleetron  phase  with  small  integration  steps  dz  =  cdt,  and  noting 
that  X  =  ct/L  so  that  dx  =  (c/L)dt  yields 


dz 

_  d^  dx  dt 
dx  dt  dz 
c  1 

-y - 

L  c 
_Vq  +  Av 
L 
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where  v  =  vo  +  Av  and  Avx4TiN(y-yo)/yo,  so  that 


f 

7  W 

Vq  -I-  471 N 

O 

1 

V 

1  ^0  )) 

L 

(4.7) 


The  simulation  is  begun  by  setting  the  electrons’  initial  conditions:  xq,  0x,  yo,  0y, 
Yo,  ^0,  and  Vq.  Using  the  equations  of  motion  (equations  (4.1)  and  (4.2)),  the 
perpendicular  positions  of  the  electrons  are  found.  From  equations  (4.3)-(4.5),  the  optical 
amplitude  and  phase  are  determined.  Next,  the  electrons’  energy  exchange  is  found  using 
equation  (4.6),  followed  by  the  electron  phase  change  from(4.7).  When  these  steps  are 
completed,  the  next  iteration  is  begun  using  the  updated  numbers,  and  the  process 
repeats. 


Below  is  an  example  output  of  a  single-mode  simulation. 


Figure  19.  Single-mode  simulation  sample  output. 

The  left  side  of  Figure  19  is  a  phase  space  plot  in  (C,  v)  similar  to  those  shown  in 
Chapter  II,  with  the  sample  electrons  shown  as  red  dots  and  the  seperatrix  shown  in  blue. 
This  plot  is  the  final  state  of  the  electrons  at  the  end  of  the  undulator,  where  the  bunching 
of  the  electrons  can  be  seen  as  well  as  the  net  decrease  in  phase  velocity.  The  upper  right 
of  the  figure  shows  the  extraction,  p,  as  the  electrons  pass  through  the  undulator  from 
T  =  0  to  1 .  The  middle  right  plot  shows  the  evolution  of  the  optical  field  strength,  a.  The 
peak  at  T  =  0.5  shows  the  nature  of  the  short  Rayleigh  length  FEL,  in  that  the  optical  field 


42 


is  focused  in  the  eenter  of  the  undulator.  The  bottom  right  plot  is  the  evolution  of  the 
optieal  field  phase  (j)  along  the  undulator. 

B.  MULTIMODE  SIMULATIONS 

Multimode  simulations  use  the  wave  and  pendulum  equations  derived  in  Chapter 
II  of  this  thesis.  An  eleetron’s  longitudinal  phase  veloeity  is  defined  by  equation  (2.38), 
whieh  again  is 

v  =  L\{k  +  kQ)P^-k\, 

where  L  is  the  length  of  the  undulator,  k  =  In/k  is  the  optieal  wave  number,  ko  =  2k/Xo  is 
the  undulator  wave  number,  and  Pz  =  Vz/  c.  For  off-axis  eleetrons,  this  phase  veloeity  is 
modified  by, 

^v  =  -[(olyl+el),  (4.8) 

where  yo  is  the  distanee  off-axis  the  eleetron  is  injeeted,  cop  =  2nNK/y  is  the  betatron 
frequeney,  and  0y  is  the  injeetion  angle,  as  in  the  previous  seetion.  [3] 

In  multimode  simulations,  the  optieal  wave  equation  is  evaluated  in  three 
dimensions  (x,y,t),  evolving  in  spaee  and  time.  The  advantage  of  the  multimode 
simulation  over  single-mode  is  that  it  ineludes  multiple  transverse  optieal  modes  and 
betatron  motion,  whieh  is  given  by 

0  _j_  0 

T  =  (To  +  ^y) cos (^  - ^/? ))  +  - -sini^co^  (r  - )) ,  (4.9) 

(Op 

where  yo  is  the  initial  off-axis  position  of  an  eleetron,  A0yo  is  the  initial  tilt.  Ay  is  a 
random  shift  due  to  the  beam  radial  spread,  A0y  is  a  random  angle  due  to  the  beam 
angular  spread,  xp  is  the  dimensionless  position  in  x  of  the  beam  foeus,  and  x  is  the 
dimensionless  time  as  deseribed  in  Chapter  II.  [2]  It  begins  with  weak  fields,  and 
evolves  over  many  passes  until  steady-state  power  is  reaehed. 

Figure  19  below  shows  an  output  from  a  multimode  short  Rayleigh  length 
simulation  with  Zq  =  0.07.  The  simulations  use  the  dimensionless  parameters  diseussed  in 
Chapter  II.  The  upper  right  seetion  lists  these  parameters.  In  this  ease,  the  dimensionless 
eurrent  density  j  =  150  and  a  dimensionless  radius  in  the  x  and  y  direetions  of  ax,y  =  0.2. 
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The  electrons  were  injected  with  an  initial  phase  velocity  Vo  =  11.25.  Both  the  electron 
beam  and  optical  mode  were  focused  at  the  center  of  the  N  =  14  period  undulator, 
denoted  by  xp  =  Xw  =  0.5.  The  beam  angular  spread  of  a9x,y  =  0.7  was  determined  by  the 
betatron  frequency  cop  =  0.6  and  the  electron  beam  radius,  and  the  Rayleigh  length  Zq  by 
the  dimensionless  mirror  radius  rm  =  2.9  and  radius  of  curvature  rc  =  0.5. 
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Figure  20.  Multimode  simulation  output. 

On  the  left-side  of  the  figure,  the  upper  and  middle  sections  show  the  optical  field 
intensity.  The  color  scale  is  defined  at  the  bottom  of  the  parameter  box,  where  dark  to 
light  blue  represents  low  to  high  amplitude,  respectively.  Sample  electrons  are  shown  in 
red  and  the  white  contour  lines  delineate  the  5%  curve  from  the  maximum  optical  field 
intensity  at  the  cavity  center.  The  upper-left  plot  shows  the  evolution  of  a  cross-section 
of  the  optical  field  amplitude  over  n  =  150  passes.  The  upper-center  plot  shows  an  end- 
on  view  of  the  optical  mode  at  the  end  of  the  last  pass. 

The  middle-left  plot  shows  the  side  view  of  the  cavity  on  the  final  pass.  The 
sample  electrons  are  shown  in  red,  and  the  optical  mode  is  shown  using  the  same  color 
scale  as  before.  Here  the  focusing  of  the  optical  field  can  be  seen,  as  well  as  the  region  of 
interaction  with  the  electrons.  The  yellow  line  at  each  end  of  the  figure  represents  the 
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profile  of  the  optieal  field  amplitude  as  it  reaehes  the  end  of  the  undulator.  In  this  ease,  in 
order  to  save  run  time,  the  simulation  shows  the  beam  only  inside  the  undulator,  from 
X  =  0  to  1 . 

The  bottom-left  plot  shows  the  evolution  of  the  eleetron  phase  veloeity  over  the 
n  =  150  passes.  The  eleetrons  are  injeeted  with  an  initial  phase  veloeity  vo  =  11.25  and 
the  phase  veloeity  spread,  and  henee  the  energy  spread,  inereases  until  the  PEL  reaehes 
saturation  at  approximately  80  passes.  The  plot  in  the  lower  eenter  shows  the  phase 
spaee  distribution  in  (^,  v)  of  the  eleetrons  at  the  end  of  the  final  pass,  demonstrating 
eleetron  bunehing. 

The  two  plots  on  the  bottom  right  show  the  evolution  of  the  wavefront  power 
(P(n))  and  gain  (ln(l+G(n)))  as  a  funetion  of  the  number  of  passes.  The  power  plot 
shows  again  that  steady  state  is  reaehed  at  approximately  80  passes.  In  this  ease, 
T)  «  3.8%  and  the  indueed  phase  veloeity  spread  was  Av  «  36.  The  extraetion  is  defined 
as 

(4,10) 

/o 

where  z  is  the  distanee  along  the  undulator.  [6] 

C.  SIMULATION  RESULTS 

Numerous  single-mode  and  multimode  simulations  were  eondueted  in  order  to 
study  the  effeet  of  ehanging  the  undulator  length,  eleetron  beam  foeus,  normalized 
Rayleigh  length,  and  output  eoupling,  and  the  results  of  the  two  simulation  teehniques 
eompared  below. 

1.  Variation  of  Undulator  Length  and  Position 

The  first  study  was  eondueted  in  order  to  determine  the  effeet  of  ehanging  the 
number  of  undulator  periods  on  extraetion.  The  undulator  wavelength  was  fixed  at 
^0  =  2.6  em  and  the  number  of  periods  varied  from  N  =  8  to  20  periods.  The  EEL  was 
below  threshold  for  N  <  8,  so  that  there  was  no  lasing.  Eigure  20(a)  below  shows  the 
extraetion  as  a  funetion  of  the  number  of  undulator  periods  for  the  single-mode 
simulation,  while  Eigure  20(e)  shows  the  results  eompiled  from  multimode  simulations. 
As  ean  be  seen  from  the  figures,  the  number  of  undulator  periods  plays  an  important  role 
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in  the  extraction  of  the  laser.  It  is  interesting  to  note,  however,  that  the  ideal  number  of 
periods  is  around  14  for  these  parameters.  Not  only  is  there  no  advantage  to  having  a 
longer  undulator  in  a  short  Rayleigh  length  PEL,  the  extraction  actually  decreases  slightly 
due  to  the  electrons  over-bunching.  When  this  happens,  electrons  actually  begin  to  take 
energy  back  from  the  optical  mode.  In  phase  space,  this  is  shown  by  the  electrons  losing 
energy  as  in  Figure  5  in  Chapter  II,  but  as  they  continue  through  the  undulator,  they 
continue  to  revolve  clockwise.  After  they  pass  the  minimum,  they  start  increasing  in 
energy,  and  are  therefore  taking  energy  back  from  the  optical  mode. 

The  difference  in  this  case  between  multimode  and  single-mode  simulations  was 
slight.  The  peak  in  extraction  in  single-mode  was  N  =  14,  while  the  multimode  yielded 
N  =  15,  both  with  a  peak  extraction  p  »  3%.  The  decrease  in  extraction  following  this 
peak  is  more  marked  in  the  single-mode  as  well. 

The  other  consideration  that  was  addressed  in  the  single-mode  simulation  is 
shown  in  Figure  20(b).  This  was  the  shifting  of  the  undulator  with  respect  to  the  electron 
beam  and  optical  mode  foci.  The  figure  shows  on  the  x-axis  where  the  electron  beam  and 
optical  mode  were  focused  in  terms  of  the  dimensionless  length  of  the  undulator.  It  is 
shown  that  the  ideal  location  for  the  foci  is  Xw=  xp=  0.5,  or  halfway  along  the  undulator. 
The  optical  mode  and  electron  beam  foci  will  be  varied  independently  in  later  sections. 
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Figure  2 1 .  Results  of  varying  undulator  periods. 

2,  Variation  of  Electron  Beam  Focus 

The  next  variation  that  was  studied  was  the  variation  of  the  amount  that  the 
eleetron  beam  was  foeused.  For  this  simulation,  the  normalized  beam  emittanee, 
determined  by  the  accelerator,  was  kept  constant  at  Sn  =  yrbOb  ~  9mm-mrad,  where  rb  is 
the  beam  waist  radius  and  0b  is  the  divergence  angle  of  the  electron  beam.  With  strong 
focusing,  then  the  radius  of  the  waist,  rb,  is  small,  but  0b  is  large.  With  weak  focusing, 
the  radius  is  large,  and  the  spreading  angle  is  small,  as  in  Figure  22  below.  This  focusing 
of  the  electron  beam  results  in  an  optimum  value. 


fj,  small,  0),  big 


rj,  big,  small 


Figure  22.  Pictorial  representation  of  electron  focusing. 

The  implication  of  this  can  be  seen  in  Figure  23  below,  which  comes  from 
multimode  output  files  similar  to  Figure  19,  where  again  the  red  is  sample  electrons  and 
the  white  line  the  5%  curve  from  the  maximum  optical  field  intensity  at  the  cavity  center. 
With  strong  focusing,  the  electrons  are  contained  entirely  within  the  optical  mode  at  the 
focus,  but  quickly  expand  outside  of  the  optical  mode,  as  in  Figure  23(a).  If  the  focusing 
is  weak,  than  the  waist  radius  will  be  larger  than  that  of  the  optical  mode,  as  in  Figure 
23(b).  This  is  important  because  for  optimal  transfer  of  energy  between  the  electrons  and 
optical  mode,  the  electrons  must  be  within  the  optical  mode. 
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Figure  23.  (a)  Strong  and  (b)  weak  electron  beam  focusing. 

Figure  24(a)  below  shows  the  extraction  as  a  function  of  the  optical  beam  waist 
radius,  rb,  for  the  single-mode  simulation.  The  extraction  peaks  because  of  the 
aforementioned  conditions  at  rb  =  0.07mm  with  a  value  of  r)  «  3.1%.  The  placement 
along  the  undulator  of  the  electron  focus  was  examined  independent  of  the  optical  beam 
as  well,  and  again  found  to  be  optimum  near  the  center  of  the  undulator,  as  shown  in 
Figure  24(b).  Figure  24(c)  shows  the  result  of  the  multimode  simulations  in  which  the 
dimensionless  beam  radius  cr  =  rf^{7t ! was  varied.  This  results  in  Omax  =  0.12 
(rb« 0.05mm)  with  an  extraction  of  p  =  3.8%. 


(a) 
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Figure  24.  Variation  of  electron  beam  waist  radius  rb. 


3,  Variation  of  the  Rayleigh  Length  zo 

In  this  section,  the  effect  of  changing  the  Rayleigh  length  will  be  examined.  The 
Rayleigh  length  was  varied  from  zo  =  1.85  cm  to  4.44  cm  with  the  cavity  length  constant, 
and  the  extraction  r|  and  mirror  intensity  I  plotted  in  Figure  25(a)  below  in  red  and  blue, 
respectively.  The  extraction  steadily  increases  as  the  Rayleigh  length  is  increased. 
However,  the  intensity  on  the  mirrors  quickly  increases  due  to  this  increased  extraction 
and  the  decreasing  spot  size.  This  means  that  the  intensity  quickly  increases  over  the 
threshold  which  the  mirrors  can  handle  before  failing.  This  level  is  assumed  to  be  200 
kW/cm  ,  and  is  shown  as  the  horizontal  green  line  on  the  figure.  Therefore,  the  Rayleigh 
length  should  be  kept  at  approximately  2.6  cm  in  order  to  avoid  mirror  damage  while 
getting  the  maximum  possible  extraction.  Figure  24(b)  again  shows  the  effect  of 
changing  the  focus  along  the  length  of  the  undulator,  this  time  moving  the  optical  mode 
waist  and  keeping  the  electron  beam  focus  at  xp  =  0.5.  The  optimum  location  for  the 
focus  is  again  shown  to  be  near  half  way  along  the  undulator. 

Figure  24(c)  shows  the  result  from  the  multimode  simulations.  In  this  case,  the 
normalized  Rayleigh  length,  which  is  normalized  to  the  length  of  the  undulator,  was 
varied.  Due  to  mirror  considerations  listed  above,  the  normalized  Rayleigh  length  should 
be  approximately  0.06,  or  2.3  cm. 
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Figure  25.  Variation  of  the  Rayleigh  length. 


4,  Variation  of  Mirror  Output  Coupling  Qn 

The  last  parameter  that  was  varied  in  both  single-mode  and  multimode 
simulations  was  the  output  coupling.  This  represents  the  amount  of  energy  that  the 
partially  transmissive  mirror  reflects  every  pass.  The  variable  in  Figure  25  below,  Qn,  is 
defined  as  the  inverse  of  the  output  coupling.  For  instance,  for  Qn  =  2,  there  is  50% 
output  coupling,  or  50%  transmission  of  the  power.  As  expected,  the  extraction  increases 
as  the  output  coupling  decreases,  which  can  be  seen  in  Figure  26(a)  and  the  blue  line  in 
(c).  This  is  due  to  the  fact  that  stronger  optical  fields  within  the  undulator  extract  more 
energy  from  the  electrons.  However,  as  was  mentioned  in  the  previous  section,  there  is  a 
limit  on  the  amount  of  intensity  that  the  mirrors  take  before  suffering  damage.  Figure 
26(b)  below  shows  the  power  contained  within  the  cavity  as  a  function  of  the  output 
coupling  for  the  single-mode  simulation,  and  is  given  by  Pm  =  QnPout-  The  pink  line  in 

Figure  25(c)  shows  the  intensity  on  the  mirrors  as  a  function  of  Qn  in  the  multimode 
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simulation.  The  green  line  in  both  plots  again  represents  the  200  kW/em  intensity  limit 
for  the  mirrors.  In  the  single-mode  simulation,  Pm  rapidly  inereases  with  Qn  and  exeeeds 
the  mirror  damage  threshold  at  a  Qn  of  approximately  2.  In  the  multimode  simulation,  it 
exeeeds  this  damage  threshold  at  a  slightly  lower  value. 


D.  MIRROR  TILT 

A  separate  set  of  multimode  simulations  was  eondueted  in  order  to  determine  the 
effeet  of  mirror  tilt  on  short  Rayleigh  length  FEL’s.  Speeifically,  the  extraetion  will  be 

shown  for  eavity  mirrors  tilted  from  0m  =  0— >0.5  with  a  normalized  Rayleigh  length  of 

0.1.  The  tilt  angle  0m  is  the  normalized  tilt  angle,  and  is  defined  as  0m  =  0  /  (/I/Tti)  . 
Other  Rayleigh  lengths  were  studied,  but  zo  =  0.1  was  the  most  instruetive  for  studying 
larger  angles  due  to  limitations  of  the  window  size  within  the  simulation.  Another 
interesting  aspect  that  will  be  addressed  is  the  effect  that  mirror  tilt  has  on  the  optimal  vo, 
which  will  be  discussed  later  in  this  section. 
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1,  Extraction  Decrease  Due  to  Beam  Tilt 

As  can  be  seen  in  the  figure  below,  the  effeet  of  tilting  the  resonator  eavity 
mirrors  is  that  the  interaction  region  between  the  optical  field  and  the  eleetron  beam  is 
not  centered  throughout  the  entire  length  of  the  undulator.  This  means  that  the  eleetrons 
are  not  seeing  the  strongest  optieal  fields  possible,  so  that  the  extraction  decreases.  The 
greater  the  mirror  tilt,  the  shorter  the  interaction  region,  and  therefore  the  weaker  the 
extraetion. 


FEI.  WAVEFRONTS 


-30  _ I, _  _ 

0  n  100  -7t/2  C,  3n/2  0  n  100 

Figure  27.  Multimode  simulation  run  with  tilted  mirrors. 


Figure  28  below  shows  the  results  of  the  runs  condueted  to  study  the  effect  of  the 
mirror  tilt  for  a  normalized  Rayleigh  length  of  zo  =  0.1.  As  ean  be  seen  in  the  figure, 
increasing  the  mirror  tilt  0m  causes  the  extraetion  p  to  deerease.  Flowever,  it  is  worth 
noting  that  the  extraetion  remains  over  1%  for  a  mirror  tilt  Bm  =  0.50.  This  corresponds 
to  an  aetual  mirror  tilt  of  0  w  0.4  prad.  With  aetive  alignment,  however,  it  is  anticipated 
that  design  toleranees  will  be  on  the  order  of  0.1  prad,  so  that  the  effect  on  extraction  due 
to  this  vibration  will  be  minimal.  [7] 
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Figure  28.  Extraction  decrease  due  to  mirror  tilt. 

2,  Mode  Shape 

Another  interesting  aspect  of  mirror  tilt  is  the  effect  on  the  optimum  value  for  vo 
and  the  mode  shape.  In  practice,  FEL’s  automatically  progress  to  the  most  efficient  value 
of  Vo  at  steady  state  for  the  given  parameters.  In  our  multimode  simulations,  that  value 
must  be  found  manually.  Iterations  of  Vo  are  conducted  and  the  maximum  extraction 
found,  which  is  what  is  plotted  in  graphs  such  as  Figure  27  above. 

For  simulations  that  do  not  involve  mirror  tilt,  as  v  is  increased  the  extraction 
increases  smoothly  to  a  maximum,  and  then  drops  sharply.  In  running  the  simulations  for 
this  study,  however,  there  was  found  to  be  dual  peaks  instead.  Both  of  these  are  shown  in 
Figure  28  below,  in  which  the  non-tilted  mirror  (0m  =  0)  is  shown  in  blue  and  the  tilted 
mirror  (0m  =  0.01)  in  red. 


Figure  29.  Example  of  finding  maximum  extraction  varying  initial  electron  phase  velocity. 
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This  is  due  to  the  ehanges  in  the  optieal  mode  shape  as  the  mirror  is  tilted.  With 
no  mirror  tilt,  the  short  Rayleigh  length  PEL  typieally  operates  in  the  fundamental 
Gaussian  mode.  With  the  mirrors  tilted,  this  mode  beeomes  distorted,  with  side  lobes 
appearing.  The  amplitude  of  these  side  lobes  increase  with  increasing  vq  until  they  reach 
a  maximum,  shown  Figure  30(a)  below,  which  shows  the  optical  mode  amplitude  profde 
at  the  cavity  center.  As  vo  continues  to  increase,  however,  the  side  lobes  are  absorbed 
into  the  main  peak,  as  can  be  seen  in  Figure  30(b).  When  this  new  mode  begins,  the 
extraction  drops.  As  vq  continues  to  increase,  the  extraction  increases  to  a  second 
maximum,  and  the  optical  mode  becomes  the  more  familiar  Gaussian,  as  in  Figure  30(c). 
This  increasing  vo  also  corresponds  to  an  increase  in  FEE  operating  wavelength. 


Figure  30.  Optical  mode  profiles  for  zq  =  0.03,  On,  =  0.01,  and  Vq  =(a)  8-0,  (b)8.5,  and  (c)10. 

In  the  case  shown  in  Figure  29,  the  FEE  evolves  from  the  “triple  hump”  mode  at 
vo»8  to  the  Gaussian  mode  at  vo«  10.5.  Knowing  Av,  the  change  in  wavelength  can  be 
found.  Equation  (2.38)  gave 


v  =  L[{k  +  k^)/3^-k]. 
Substituting  the  definition  k  =  2n/X  yields 

v(X)='^Lp,*KLIi.~L. 

Taking  the  change  in  v  as  a  function  of  X  yields 

Av_  XtvP^L  ItvL 


(4.11) 


(4.12) 
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Equation  (2.28)  showed  that 


and  equation  (2.16) 


Combining  these  two  yields 


Inserting  equation  (3.19)  into  (3.18)  yields 


A  V  2nL  f .  .  /I  ^ 
- =  — ^  1-1  +  — 

AA  ^  V  ^0  y 


(4.13) 


(4.14) 


(4.15) 


so  that 


IttL  AA 
Av  = - . 

Aq  a 


Sinee  L  =  NXq,  then  equation  (3.21)  ean  be  written  as 

AA  1  , 

—  = - Av  . 

A  InN 


(4.16) 


(4.17) 


For  A^=  18  and  Av  ~2.5,  equation  (3.22)  shows  that  the  wavelength  ehange  in  the  FEE 
due  to  this  transition  is  approximately  AAIX  «  2.2%. 
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V.  CONCLUSION 


The  free  eleetron  laser  is  a  promising  new  weapon  system  that  will  ehange  the 
faee  of  modern  warfare.  Speed  of  light  delivery  of  energy  on  target,  deep  magazines,  low 
eost-per-shot,  and  pinpoint  aeeuraey  offer  signifieant  advantages  over  eurrent  weapon 
teehnologies.  Although  the  eurrent  eoneentration  for  PEL  applieations  is  on  ship  self 
defense,  there  are  numerous  other  applieations  mentioned  in  the  introduetion  of  this 
thesis,  and  undoubtedly  others  that  will  be  realized  onee  it  is  implemented  in  an 
operational  environment. 

Simulations  eondueted  by  the  Direeted  Energy  and  Eleetrie  Weapons  Center  at 
the  Naval  Postgraduate  Sehool  show  that  a  short  Rayleigh  length  is  an  important 
eonfiguration  for  a  eompact,  weapons  grade  EEE  to  reduee  intensity  on  mirrors  in  the 
resonator  eavity.  Due  to  the  limited  spaee  allowed  for  a  weapon  system,  the  shorter 
resonator  eavity  lends  itself  to  ship-board  implementation.  The  simulations  eondueted 
for  the  2003  International  Tree  Eleetron  Easer  Conferenee  in  Tsukuba,  Japan,  whieh  have 
been  ineorporated  into  this  thesis,  define  some  of  the  parameters  neeessary  for  next 
generation  EEEs.  Optimal  values  for  parameters  sueh  as  the  number  of  undulator 
periods,  eleetron  beam  foeus,  the  normalized  Rayleigh  length,  and  mirror  output  eoupling 
were  determined.  Simulations  were  also  eondueted  in  order  to  examine  the  effeet  of 
mirror  tilt  on  laser  power  and  extraetion.  This  demonstrated  the  importanee  of  aecurate 
simulations,  in  that  the  effeet  of  varying  parameters  ean  be  determined  in  order  to  help 
seientists  and  engineers  develop  future  systems,  and  eomparison  between  the  single¬ 
mode  and  multimode  simulation  teehniques  eontributed  to  validating  the  results. 
Simulations  also  aid  greatly  in  understanding  how  EEE  parameters  interaet. 

Multiple  beam  direetors,  already  proposed  for  maximum  radial  eoverage,  ean 
mitigate  thermal  blooming.  The  power  limitations  eaused  by  thermal  blooming  ean  be  a 
signifieant  faetor  in  the  implementation  of  high-power  lasers.  The  use  of  multiple 
direetors  will  allow  for  more  power  to  travel  through  the  atmosphere  beeause  eaeh 
individual  beam  will  be  below  the  threshold  of  thermal  blooming,  while  the  eombined 
intensity  on  target  still  meets  the  desired  speeifieations.  The  study  eondueted  in  this 
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thesis  quantifies  the  effeets  of  beam  direetor  separation  on  thermal  blooming,  and  shows 
that  this  is  a  viable  option  for  the  reduetion  of  thermal  blooming. 

Direeted  energy  weapons  promise  to  revolutionize  modem  warfare.  The  results 
of  this  thesis  show  how  a  free  eleetron  laser  ean  be  designed  and  deployed  on  a  ship  for 
use  as  a  defensive  weapon,  addressing  key  issues  sueh  as  mirror  damage,  eavity  stability, 
and  thermal  blooming. 
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